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We use the lattice model of directed walks to investigate the conforma- 
tional as well as the adsorption properties of a semiflexible homopolymer chain 
immersed in a good solvent in two and three dimensions. To account for the stiff- 
ness in the chain we have introduced energy barrier for each bend in the walk 
and have calculated the persistent length as a function of this energy. For the 
adsorption on an impenetrable surface perpendicular to the preferred direction 
of the walk we have solved the model exactly and have found the critical value 
of the surface attractions for the adsorption in both two and three dimensions. 
We have also enumerated all the possible walks on square and cubic lattices for 
the number of steps TV < 30 for two-dimensions and N < 20 for three dimen- 
sions and have used ratio method for extrapolation. The transition located using 
this method is in excellent agreement with the results found from the analytical 
method. 

PACS numbers: 64.60.-i,68.35.Rh,05.50.+q 

I. INTRODUCTION 

Biopolymers are known to exhibit under different environments a variety of persistent 
lengths ranging from being much smaller than the over all length of the polymer, to being 
comparable to the chain length |]J . When the persistent length associated with the polymer 
is much smaller than the overall length of the chain, the polymer is said to be flexible. On 
the other hand, when the persistent length is comparable to the chain length, the polymer 
is said to be rigid. When the persistent length falls in between the two extremes, the chain 
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is said to be semiflexible. The conformational properties of such chains have attracted 
considerable attention is recent years because of experimental developments in which it has 
become possible to pull and stretch single molecule to measure elastic properties 0. Such 
studies reveal a wealth of informations about the conformational behaviour of semiflexible 
polymers that are of clear biological importance. 

An impenetrable surface is known to affect the conformational properties of polymers in a 
significant way H |J. This is due to a subtle competition between the gain of internal energy 
and a corresponding loss of entropy at the surface. Since the flexibility of the chain affects 
this competition, a semiflexible chain is expected to show different adsorption behaviour 
compared to that of a flexible polymer chain. A stiff chain is known to get adsorbed easily 
compared to a flexible chain ||. 

A simple way to account the stiffness of a semiflexible chain is to constrain the angle 
between the successive segments to be fixed. The value of the angle depends on the local 
stiffness of the chain. This prescription leads to the freely rotating chain model ||. In 
the continuum limit the freely rotating chain becomes the so called worm like chain (WLC) 
0. In these models the persistent length l p is defined as a characteristic length for tangent- 
tangent correlation function < t(s)t(s) >~ exp(-^j^). The tangent vector t(s) is defined 
as , where r(s) is parametrized in terms of the arc length s of the chain 0. 

Though the worm like chain model of Kratky and Porod has been used extensively 
to study the conformational properties and surface adsorption of a semiflexible chain ||, it 
can not mimic exactly the dimensional behaviour of the real chains. In this paper we use 
the lattice model of directed walk and introduce stiffness in the chain by associating energy 
with every bend of the walk and calculate the bulk and adsorption properties of the chain 
function of stiffness of the chain. 

The paper is organized as follows: In Sec. II we describe the lattice model of directed 
walk and investigate the bulk properties. We calculate the value of persistent length as a 
function of energy associated with the bend. In Sec. Ill we discuss the surface adsorption of 
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a semiflexible chain represented by a directed walk on a plane perpendicular to the preferred 
direction of the walk. We also use the exact enumeration technique to locate the adsorption 
desorption transition and compare the result with those found exactly. This allows us to 
comment upon the accuracy of the method of exact enumerations. 

II. A LATTICE MODEL FOR THE SEMIFLEXIBLE CHAIN 

We consider a model of a self-avoiding directed walk on a lattice [[|. Though the 
directedness of a walk amounts to some degree of stiffness as all directions of the space is 
not treated equally, stiffness in the chain is introduced by associating energy barrier with 
every turn of the walk. Though the model is very restrictive in the sense that the bend can 
be either 90° or no bend at all, The model can be solved analytically and therefore gives 
the exact values of conformational and adsorption properties of a semiflexible chain. We 
consider two specific cases of directed walks: If the walker is allowed to take steps along 
±2/-axis (in two-dimensions) and only along +x-axis the walk is said to be partially directed- 
self-avoiding walk (PDSAW). On the other hand, if the walker is allowed only along +y 
and +x directions then the walk is said to be fully directed-self-avoiding walk (FDSAW). 
In the case of three dimensions {3D), a PDSAW is one in which walker is allowed along in- 
direction but only along +x and +z directions while in the FDSAW the walker is allowed 
to move only along +x, +y and +z directions. 

A stiffness weight k = exp(—f3eb) where (3 = (fe^T) -1 is inverse of the temperature and 
e&(> 0) is the energy associated with each turn. For k = 1 or e& = the chain is said to be 
flexible and for 0<A;<lorO<eb<oo the chain is said to be semiflexible. When % — > oo 
or k — > 0, the chain becomes a rigid rod. 

The partition function of such a chain can be written as 

=£*ir e * N k Nb (i) 

allwalksof N steps 

Here N b is the total number of bends in a walk of N steps and x is the step fugacity. 

3 



A. Conformational properties of the chain in Two-Dimensions 



(i)The case of PDSAW: 

In 2D, the generating function for PDSAW has two components; one along the direct- 
edness {i.e. + x-axis) and other perpendicular to it (i.e. ± y-axis) as shown in Fig.l. The 
recursion relation for these two components of generating function are [[J; 

X = x + x(X + 2kY) (2) 



Y = x + x(kX + Y) 



(3) 




Fig. 1; The diagrammatic representation of the recursion relations Eqs. (2) and (3), for 
PDSAW . The thick arrows X and Y denote all possible walks with the initial step along +x and 
±y directions respectively. 



Solving Eqs. (2) and (3) we get 



X 



Y 



x+(2k-l)x 2 
l-2x + x 2 - 2x 2 k 2 



x + (k — l)x 2 



1 — 2x + x 2 — 2x 2 k 2 
The partition function can therefore be written as 



(4) 



(5) 



>™ d ( X}k )=X + 2Y=- [ ± 



(4Jfe - 3)x 2 + 3x 



2x + x 2 - 2x 2 k 2 

The critical point for polymerization of an infinite chain is found form the relation 



(6) 



1 - 2x + x z - 2x 2 k 2 = 



(7) 



This leads to the critical value of the step fugacity for a given value of k as x c = 



1+V2k' 




ID fully directed 
2D Partially directed 
3D partially directed 
3D fully directed 



0.2 



Fig. 2; The variation of step fugacity x c with 0e b . 

The stiffness in the chain increases the value of fugacity for polymerization. This de- 
pendence is shown in Fig. (2) by long-dashed line in which we plot function of (5e b . 
We define the persistent length as the average distance between two successive bends of the 
walk, i.e. 



l p =< L> I < N b > 



(8) 



Where L =< N > a, a being the lattice parameter. 
For the PDSAW in 2D we find 

l P =^^[V2 + exp((3e b )\ 

The dependency of l p on (3e b is shown in Fig. (3) by a long-dashed line. 



(9) 
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The value of l p increases exponentially with the bend energy at a given temperature, 
(ii) The case of FDSAW: 

In this case the polymer is directed along +x and +y direction; leading to the following 
recursion relations for the generating functions, 

X = x + x(X + kY) (10) 



Y = x + x(kX + Y) (11) 

Solving these equations we get the following value for the partition function Zj d d ; 

Zf d =X + Y = 2X - (12) 

M l-(l + k)x K J 

The critical value of step fugacity is found to be x c = The variation of x c with ftey, 
is shown in Fig. (2) by solid line. 

The value of persistent length in this case attains a simple relation, i.e. 

l p = l + e^ b (13) 

The variation of l p with the bending energy is shown in Fig. (3) by a solid line. 
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B. Conformational properties of the chain in Three-Dimensions 



(i)The case of PDSAW: 

In the case of PDSAW the polymer chain is directed in two-directions. The recursion 
relations for generating functions are || 

X = x + x(X + 2kY + kZ) (14) 
Y = x + x(kX + Y + kZ) (15) 

Z = x + x(kX + 2kY + Z) (16) 

Solving these equations we get the values of X, Y, Z and the partition function as 

Y _v_ x+(2k-l)x 2 

(1 + k - Ak 2 )x 2 - (k + 2)x + 1 1 ; 

Y = * + {k-l)x 2 (18) 

(1 + k - Ak 2 )x 2 - (k + 2)x + 1 v ; 

Z pA _ -X+2Y + Z- {1 + k _ 4k2)x2 _ {k + 2)x+1 (19) 

That is at x c = 2{k+i^H^) ' ^ e ^p.d diverge. In this case the dependence of the 
fugacity for polymerization on the stiffness is more involved compared to the case in 2D. 
The variation of x c with /3e b is shown in Fig. (2) by the dashed line. For persistent length 
we find 

_ [85 + 19^17 - (102 + 26^/17) exp((3e b ) + (34 + S^Y!)exp{2f5e b )\2[exp{2(3e b ) + exp((3e b ) - 4] 

Lp 



(1 - V17 + 2exp{(3e b ))[2$A + 52V17 - (272 + QAV17)exp((3e b )] + (85 + 2Wl7)exp(2/3e b ) 

(20) 

The value of l p as a function of f3e b is plotted in Fig. (3) by the dashed line. 



(ii)The case of FDSAW 

In this case the polymer is directed along all the three directions i. e. along +x, +y and 
+z directions. We can write following recursion relations 

X = x + x(X + kY + kZ) (21) 

Y = x + x(kX + Y + kZ) (22) 

Z = x + x(kX + kY + Z) (23) 

The solution of these equations leads to 

X = Y = Z = , X - (24) 

Thus the partition function of the system can be written as 

Zr i =X + Y + Z= T - JTTWx (25) 

The critical value of the step fugacity is x c = (Y+2fcj ■ The variation of step fugacity with 
bending energy (3e b is shown in Fig. (2) by a dot-dashed line. In this case l p is found to be 

l P = 1 + \exp{Pt b ) (26) 

The value of l p as a function of (3e b is plotted in Fig. (3) by a dot-dashed line. 

In all the cases discussed above the persistent length shows exponential dependence on 
the bending energy. 

III. SURFACE ADSORPTIONS 

In the case of directed model we have two distinct surfaces; one parallel and the other 
perpendicular to the directedness of the walk. The adsorption of polymer on a surface 
parallel to the preferred direction of the walk has been studied in case of 2D using the 
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transfer matrix method 0. The features associated with the adsorption were found to be 
same as in the isotropic case except that the critical value of surface attraction for adsorption 
is higher. The surface perpendicular to the direction of walk may give different features as 
the walk once leaves the surface it can not return to it due to restriction on the walk. Here 
we report the result found analytically for the adsorption of directed semiflexible chain on 
a surface perpendicular to the directedness of the chain both in two and three dimensions. 

A. Adsorption of a directed semiflexible chain on a surface perpendicular to the 



(i)The case of PDSAW: 

In case of two- dimensions, surface is a line represented by x = 0. Let S be the component 
of generating function along the surface and X the component perpendicular to the surface 
as shown in Fig. (4). Following the method outlined above we can write surface component 
as 



where s = ujx, and uj = exp(—/3e s ) being the weight associated with each step along the 
wall. 



Fig. 4; The diagrammatic representation of the recursion relation (28). Each walk of the 
polymer chain starts from O. In this digram X and S denotes all possible walks with initial step 
along the +x and along the wall respectively. 



directedness of the chain in two-dimensions 




(27) 




9 



For uj = 1 Eq. (27) reduces to Eq. (5). 

The partition function in presence of surface for PDSAW is found to be 

Z™ lL (k,w,x)=X + 2S 

Combining with Eq. (4) we find 

2d 2sx 2 (l - 2k) + 2s{l - 2x) + (2sk + 1 - s)[x + (2k - l)x 2 } 
SpA - ~ (I - s)(l - 2x + x 2 - 2k 2 x 2 ) 

The critical value of adsorption transition is found from the relation 



1- s)(l-2x + x 2 -2k 2 x 2 ' 







(28) 



(29) 



(30) 



This leads to u c = ^- = \[2k + 1, which reduces to lo c = \/2 + 1 for the flexible polymer 
chain. The variation of u c with /3e& is shown in Fig. (5) by a long dashed line. 



2D fully directed 
2D partially directed 
3D partially directed 
3D fully directed 




Fig. 5; The exact value of iv c for different values of f3eb- The lines in this figure correspond 
to analytical results however the dots on the lines correspond to the value obtained from exact- 
enumeration method in 2D and the cross used to denote the uj c value for 3D partially directed 



case. 



(ii)The case of FDSAW: 

It is straight forward to show that for FDSAW the partition function in the presence 
of surface is 
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S = s + s(s + kX) + s 2 (s + kX) + s 3 (s + JfeA") + . . . (31) 



Using value of X from Eq. (12) we get 

zr = x + s= s(l-(l + k).) + .( S k-s +lj 

(1 _ s )(l _ (l + k)x) 

Which gives adsorption transition point at uj c = 1 + k. The variation of uj c with fie^ is 
shown in Fig. (5) by a solid line. 

B. Adsorption of a semiflexible directed chain on a surface perpendicular to one out 
of the two preferred direction of the chain in three dimensions 

(i)The case of PDSAW: 

The analysis given above can be generalized in 3 dimensions where surface dimension is 
two i.e. x — y plane at z — 0. In the case of PDSAW as mentioned above the choice of the 
walker is restricted to the +x-axis, ±y-axis and +z-axis. Let S x and S y is the component 
of the total partition function Z^ d along +x and ±y axis respectively, however component 
perpendicular to the wall along +z axis remains same as defined by Eq. (17). We can, 
therefore write 

_ s - s 2 + 2s 2 k + Z(2k 2 s 2 + sk- s 2 k) 
bx ~ 1 - 2s + S 2 _ 2s 2p ( s < V ^ 

_ s - s 2 + s 2 k + Z(k 2 s 2 + sk- s 2 k) 
y ~ 1 - 2s + s 2 - 2s 2 P { } 

The expression for the partition function in this case found from the relation 

Zsl d d . = S x + 2S y + Z (35) 

Substituting the value of S x , S y and Z we find 

M x (l- x + 2kx)U + S (3-3s + 4sk)V 

SpAX ' ' ] (l-2s + s 2 -2s 2 A; 2 )[(l + A;-4A; 2 )a; 2 -(A; + 2)x+l] 1 ] 

Where U and V are, 
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U = 1 - 2s + 3sk + s 2 - 3s 2 k + 2s 2 k 2 



(37) 



V = 1 - 2x - kx + x 2 + kx - 4k 2 x 2 (38) 



The two singularities appearing in eqn.(37) give the critical value of x c = 2{i+k-Ak 2 anc ^ 
uj c = n^^rf^ViTfc) ' ^ or fl ex ibl e polymer chain (i.e. k = 1 or e& = 0) it gives x c = ~ 3+ ^^ 
and u c = 1.47524 

(ii)The case of FDSAW: 

For FDSAW, the partition function can be easily be evaluated. Here we write the final 
form of the partition function as 

Z, % .( k ,», ,) =S . + 5, + Z = (1 l * + _'-fZ - (2 l S f 1)x) (39) 
The two singularities appearing in Eq. (39) gives the critical value of x c = ^ttt and 



2k+l 

2fe+l 



c fe+1 • 

The variation of uo c with /3e& are given in Fig. (5) for PDSAW and FDSAW by a dashed 
and dot-dashed lines respectively. 

IV. RESULT FROM EXACT ENUMERATION METHOD 

Since the analytical approach is limited to very few cases, one often has to resort to 
numerical methods, such as Monte Carlo simulations or a lattice model using extrapolation 
of exact series expansion (referred to as exact enumeration method). The later method has 
been found to give satisfactory results as it takes into account the corrections to scaling. 
To achieve the same accuracy by the Monte Carlo method, a chain of about two orders of 



magnitude larger than in the exact enumeration method has to be considered flO|. 

We have enumerated all possible walks of length N < 30 on square lattice and of length 
N < 20 on cubic lattice. The canonical partition function is written as 

Z N (k, ^) = E E C n(N s , N b )iu N °k N >> (40) 

N s N b 
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Here N s is the number of steps on the surfaces. The reduced free energy per monomer 
is found from the relation 

G(k, u) — lim —logZj^(k,uj) (41) 



The limit iV — > oo is found by using the ratio method |TJJ for extrapolation. 



The transition point for adsorption-desorption is found from the maximum of 9 (= 
9< fe s> )- The transition points found from this method are shown in Fig. (5) by dots and 
cross. The results found from this method are in very good agreement with those found 
exactly in above sections. This result indicates that as for as locating the adsorption- 
desorption transition of a long flexible as well as semiflexible chains immersed in a good 
solvents are concerned the method of exact enumeration can give reliable results. 



V. CONCLUSION 

In spite of the sever restriction imposed on the angle of bending of the chain, the lat- 
tice models may provide interesting results for the conformational and surface adsorption 
properties of a semiflexible chain. Introducing directedness in the walk allowed us to solve 
the model exactly in both two and three dimensions. We have calculated the step fugacity 
for polymerization of an infinite chain and the persistent length as a function of bending 
energy associated with bending. We have also been able to obtain to the critical value for 
adsorption of a directed chain on a surface perpendicular to the preferred direction of the 
walk analytically in both two and three dimensions. The dependence of this critical value 
of surface attraction on the stiffness of the chain have been evaluated. 

We have also examined the accuracy of the method of exact enumeration in locating the 
adsorption-desorption transition and have found that the method give values that are in 
excellent agreement with the exact values. 
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